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ABSTRACT 

We argue that a universal density profile for dark matter halos arises as a natural 
consequence of hierarchical structure formation. It is a fixed point in the process of 
repeated mergers. We present analytic and numerical arguments for the emergence 
of a particular form of the profile. At small radii, the density should vary as r _Q , 
with a determined by the way the characteristic density of halos scales with their 
mass. If small halos are dense, then a is large. The mass-density relation can be 
related to the power spectrum of initial fluctuations, P{k), through 'formation time' 
arguments. Early structure formation leads to steep cusps. For P(k) ~ k n we find 
a ~ 3(3 + n)/(5 + n). The universal profile is generated by tidal stripping of small 
halos as they merge with larger objects. 

Key words: cosmology: theory - dark matter - stellar dynamics - galaxies: kinemat- 
ics and dynamics 
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1 INTRODUCTION 

In a series of recent papers, Navarro, Frenk & White (1995, 
1996, 1997, hereafter NFW) have claimed that dark halos 
formed by dissipationless hierarchical clustering from gaus- 
sian initial conditions can be fitted by a universal density 
profile of the form 



p(r) = 5- 



R 



(R + : 



(1) 



where R is a characteristic length scale and 8 a characteris- 
tic density. This function fits the numerical data presented 
by NFW over a radius range of about two orders of mag- 
nitude. Equally good fits are obtained for high mass (rich 
galaxy cluster) and for low mass (dwarf galaxy) halos, and 
in cosmological models with a wide range of initial power 
spectra P(k), density parameters Q, and cosmological con- 
stants A. In any given cosmology the simulated halos show 
a strong correlation between R and 5; low mass halos are 
denser than high mass halos. NFW interpret this as reflect- 
ing the fact that low mass halos typically form earlier. In 
independent work Cole & Lacey (1996) came to similar con- 
clusions based on studies of a series of cosmological simula- 
tions with P oc k" and £7 = 1. We shall refer to the profile 
of equation (1) as the NFW profile. 

It will prove useful below to consider a broader family 
of density profiles: we will adopt the family defined by 



p{r) 



(2) 



This family has been studied extensively by Zhao (1996); 
subfamilies were studied by Dehnen (1995) and by Tremaine 
etal. (1994). The NFW profile is of this form: an inner cusp 



with logarithmic slope a — 1, an outer envelope with loga- 
rithmic slope j3 — 3, and a 'turn over exponent' 7 = 1. An- 
other well-known example is the Hernquist (1990) profile: 
(a,/3,7) = (1,4,1). Both NFW and Cole & Lacey (1996) 
compared the Hernquist profile to their numerical data. Al- 
though it fits high mass halos almost as well as the NFW 
profile, its more rapid fall-off does not fit low mass halos 
adequately. Detailed comparisons of both profiles with high 
resolution simulations of galaxy clusters are also described 
by Tormen, Bouchet & White (1996). 

In this paper we argue that a universal profile arises 
as a result of repeated mergers. Violent relaxation of a fi- 
nite mass isolated system leads to a profile with p oc r -4 
as r tends to infinity (Aguilar & White 1986, Jaffe 1987, 
Merritt, Tremaine & Johnstone 1989, Barnes & Hernquist 
1991). Sufficiently strong violent relaxation, either through 
repeated merging or through cold inhomogeneous collapse, 
produces p oc r -3 over the regions which contain most of 
the mass (White 1979, Villumsen 1982, Duncan, Farouki & 
Shapiro 1983, McGlynn 1984). There are, however, good rea- 
sons why neither of these results should apply to halo forma- 
tion through hierarchical clustering. Cosmological mergers 
are not isolated, they occur continually and from bound or- 
bits, and they usually involve objects of very different mass. 
These properties are reflected in the other simple paradigm 
for halo formation, the spherical infall model (Gunn & Gott 
1972, Fillmore & Goldreich 1984, Hoffman & Shaham 1985, 
White & Zaritsky 1992). This model predicts relatively shal- 
low density profiles, p ~ r~ 2 for the virialised regions of 
halos. 

The next section uses simple arguments to show why 
merging in a cosmological context might give rise to a char- 
acteristic density profile with a central cusp. The third sec- 
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tion then carries out some numerical experiments to test 
these arguments. A final section discusses our results. 



2 MERGERS AND TIDAL DISRUPTION 

Hierarchical clustering proceeds through a continual series 
of unequal mergers. Small objects typically form at earlier 
times and so have higher characteristic densities. The results 
of Navarro, Frenk & White (1996) for a CDM universe can 
be fit by a power law 



S oc M 



(3) 



with v ~ 0.33. The characteristic radius is given by 

i?«M (1+ " )/3 . (4) 

Simple scaling arguments using linear theory (Kaiser 1986) 
predict that if P(k) ~ k" then 



v= (3 + n)/2. 



(5) 



v — 0.33 corresponds to n = —2.3 as expected for a CDM 
universe on the relevant scales. 

Consider a large parent halo merging with a smaller 
satellite system. Two important dynamical processes are re- 
sponsible for the evolution of the system. Dynamical friction 
causes the satellite's orbit to decay, bringing it down into 
the central regions. Tidal stripping removes material from 
the satellite and adds it to the diffuse mass of the parent. 
Our claim is that the combination of these processes leads 
to a ~ 1 in the inner cusp. This structure is a stable fixed 
point in the process of repeated mergers, in a sense to be 
defined below. 

A crude but useful way of understanding tidal strip- 
ping, going back at least to the work of van Horner (1957), 
is as follows. Consider a satellite of mass m orbiting a par- 
ent of mass M at mean distance r. The frequency f2 of the 
orbit is given roughly by Q, 2 ~ Pm{t), where Pm{t) is the 
average density of the parent interior to r. Next consider a 
particle Q orbiting within the satellite at a mean distance 
£ from its centre. The frequency uj of this orbit is given by 
~ Pm(£), where Pm(£) is the average density of the satel- 
lite interior to £. If Q = lo there is a resonance between the 
force the satellite exerts on Q, and the tidal force exerted 
by the parent. This results in a secular transfer of energy 
to Q, which then escapes from the satellite. Thus, crudely 
speaking, the satellite is stripped of material down to ra- 
dius £ defined by w(£) = Q(r) . The newly stripped material 
will orbit throughout the parent but will, on average, re- 
main near the orbit of the satellite at the instant that it was 
stripped. 

Weinberg (1994abc, 1996) has shown that material can 
actually be stripped off the satellite substantially within the 
primary resonance Q. = u>. The strongest effect is neverthe- 
less at this resonance, and since the timescale of real cosmo- 
logical mergers is not very different from there is little 
time for weaker effects to be felt. For simplicity we stick to 
the standard assumption of stripping at the primary reso- 
nance, and we will use iV-body experiments to check our 
principal results. 

Since clustering is hierarchical we expect the satellite to 
have a higher characteristic density and a smaller character- 
istic radius than the parent. As a result the density of the 



a < 3jV(1 + v) 




Rm 



Figure 1. The superimposed density profiles of a parent halo of 
mass M and of homologous satellite halos of mass m < M. The 
logarithmic slope of the central cusp is —a. Two satellites are 
shown, one with p m {Rm) < dm (lower solid curve) and one with 
pm(Rm) > <Sm (upper solid curve). The former has 3v/(l+v) < a 
and the latter 3u/(l + v) > a. 

parent at the characteristic radius of the satellite pM(Rm) 
can be either larger or smaller than the density of the satel- 
lite at the same radius p m (Rm) = S m (see Figure 1). (For 
a < 3 the average density is also proportional to r~ a .) Thus, 
when dynamical friction has brought the satellite right to 
the centre of the parent, it will either have been entirely 
disrupted, or will have survived largely intact; the result de- 
pends on whether pM(Rm) is less than or greater than 5 m - 
The case of equality, where the satellite is marginally dis- 
rupted by the parent as it reaches the centre, is given by 



Pm(Rti 



(6) 



which can be written in terms of the masses using equations 
(4) and (3) as follows 



1=(£) 

whence 
3^ 



a ■ 



1 + v 



In conjunction with (5) this predicts 
\5 + nJ 



(J) 



(8) 



(9) 



It is interesting that this exponent is the same as that pre- 
dicted in a completely different context, but for similar rea- 
sons, for the small scale behaviour of mass autocorrelation 
function in hierarchical clustering on the stable clustering 
hypothesis (Davis & Peebles 1977). These equations consti- 
tute our prediction for the logarithmic slope of the inner 
cusp in cosmological halos. For the Navarro, Frenk & White 
(1996) halos, v ~ 0.33 implying a ~ 0.75. This is probably 
consistent with NFWs fits to these halos since their exper- 
iments do not resolve the central cusp to r <C R- In fact 
a does not vary strongly with n for cosmologically interest- 
ing values of —1 to —2.5. The simulations of Cole & Lacey 
(1996) include the cases n = {0, —1, —2} corresponding to 
a = {1.8, 1.5, 1}. They do not resolve the central cusps well, 
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but there is a hint in their Figure 9 that n — has steeper 
halos than n — — 2. 

A further argument is needed to show that halos will 
actually prefer the special value a = 3v/(l + v). Suppose a 
in a parent halo were smaller than this value, and it merged 
with a homologous satellite. In this case pM^Rm) < S m , 
so the satellite survives intact right to the centre. It thus 
boosts the density of the parent cusp in the region r < R m , 
and steepens the effective logarithmic slope in r < Rm', 
i.e. a increases. Now suppose a in the parent were larger 
than 3u/(l + v). In this case pM{Rm) > 8 m , so the satellite 
will be disrupted at r > R m , spreading its mass out over a 
region which is larger than its own characteristic radius. It 
thus boosts the density of the cusp in the region r > R m , 
and softens the effective logarithmic slope in r < Rm ■ Thus 
hierarchical halo formation by repeated mergers should lead 
to the density profile taking on a fixed point form, which 
depends on the mass-density relation of merging subunits. 



1.5 



0.5 








5 10 15 

iteration 

Figure 2. The logarithmic slope of the cusp a as a function of 
the number of mergers in the experiments of Section 3.1. From 
bottom to top at the right v = (0.33, 0.5, 0.75). 



NUMERICAL EXPERIMENTS 



We now test the assertion that there is a fixed point in the 
process of repeated mergers. In outline our procedure is: 
to merge a halo repeatedly with a small satellite that is 
created by scaling the parent according to equation (3). We 
emphasise that the satellite is always a copy of the current 
halo, i.e. the latest merger product. 

First we take a parent halo with mass M but with an 
initial density profile quite different from the NFW profile: 

R? 

P^= \ R 2 +r 2y/2 > ( 10 ) 

with 5 = 1, R = l and /3 an input parameter. The following 

steps are then carried out repeatedly: 

1 Choose a mass ratio for the next merger, / = m/M, 
by taking a random value from a distribution between 
/ min = iQ- : i/( 1 +^ an d / max = 0.4 in which the number 
of daughters with mass ratios in (/, / + d/) is specified 
by 



dN = pdf. 



(11) 



Thus 7] — gives a uniform distribution of daughter 
masses, and r\ < produces a larger number of low 
mass daughters. 

Create a satellite of mass fM by cloning the parent 
homologously: p -*■ pf~ v , r — » r / (1+ " )/3 . 
Merge the satellite with the parent. Section 3.1 de- 
scribes a semi-analytic model of the merging process 
based on the arguments in Section 2. Section 3.2 de- 
scribes an JV-body implementation of the merging pro- 



3.1 Semi-analytic models 

First we implement our simplified picture of tidal stripping 
as described in Section 2. On a logarithmic grid in r between 
r — 0.01 and r = 100, a merger is simulated as follows. Work 
inwards from the outer edge of the parent, removing all the 
mass from the satellite where p m (£) = Pm(t) and placing it 
in the radial bin at r. If the profiles cross (pM(R m ) < dm) 
add the satellite bin by bin to the parent, as if it survived 
intact to the centre. 




P0.001 



0.00001 
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Figure 3. The parent halo profile from the last iteration in the 
experiment with (y, rj, 13) = (0.33, 0, 4) superimposed on a homol- 
ogous satellite with mass ratio / = 0.22. Compare Figure 1. 



Figure 2 shows the evolution of the profile as the iterations 
proceed for runs with (/?, (3) = (0, 4) and v = (0.33, 0.5, 0.75). 
Clearly the homogeneous core is steepened and settles down 
into a cusp, as it was argued it should in the previous section. 
A steady state is reached, which, as predicted, automatically 
flattens itself if it gets too steep. The final value of a, av- 
eraged over 10 iterations, is summarised in Table 1 and in 
Figure 4. As shown by Figure 4, the value of a is systemat- 
ically larger than predicted by equation (8). Figure 4 shows 
the results of experiments with (3 = 3 and j3 = 4. The lower 
value of f3 tends to produce slightly lower values of a, but 
still consistent with the f3 = 4 results at the la level. 

Figure 3 shows the profile of a satellite superimposed 
on its parent, so it is the equivalent of Figure 1. Since a is 
slightly greater than predicted by equation (8), the satellite 
is less dense than its parent at small radii. Note that the 
inner cusp is very close to a power law. However, Figure 5 
shows that the slope a is correlated with the turn over expo- 
nent 7. This could mean that we are not resolving the power 
law in the cusp. To test this we compare the measured value 
of a with the actual logarithmic slope of the profile at the 
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Table 1. The results of the experiments of Section 3.1. The input 
parameters are (/3,?7, v), and the output fitted profiles are spec- 
ified by a and 7. The errors in the last two decimal places are 
given in brackets (la estimated from the last 10 iterations in each 
case) . a e g/a is a measure of how well we resolve the inner cusp 
(see text). 
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0.62(08) 
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0.7 


3 





0.2 
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0.75 


1.47(14) 
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Figure 4. The relationship between v and a as measured (with 
la error bars). Times signs are for (3 = 4, squares are for (3 = 3 
(both with r) = 0) and no symbol for (i = 4 and different values 
of r/ (see Table 1). The solid line is the prediction of equation (8). 

inner grid point: 

_ dlnp a + Pr^'< 
Qcff -"dln7 = T+r^h ■ (12) 

The last column of Table 1 shows that the difference between 
a and a e g is small, thus we have adequately resolved the 
central power law. 

Calculations were performed with ((3, v) = (4, 0.33) and 
/ drawn from different distributions. We adjust the relative 
numbers of large and small mergers by changing the value of 
77 in equation (11). The results are summarised in Table 1, 
and also appear in Figures 5 and 4. Generally the result- 
ing values of a are similar, show no trend with rj, and are 
statistically consistent with each other. There is nothing in 
the argument leading to equation (8) which depends on the 
mass ratio, so this is expected. The measured 7 also shows 
no distinct trend with rj, although, as can be seen from Fig- 
ure 5, varying r\ can lead to large scatter in the value of 7. 
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Figure 5. The correlation between 
Symbols as in Figure 4. 
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7 and a (with la error bars). 
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Figure 6. From bottom to top the ratio of the 10 — 50 and 
80, 70, 60%-mass radii to the half-mass radius of the halo in sim- 
ulation A as a function of the number of mergers. Each quantity 
is normalised to its initial value. 

All these values are, however, consistent with each other at 
the la level. 



3.2 A-body experiments 

The code we used was the publicly available adaptive P 3 M 
code of Couchman, Pearce and Thomas (1996) (configured 
to simulate isolated, collisionless systems). Compared with 
these A^-body experiments, the simulations of Section 3.1 
have better resolution and they are much cheaper to per- 
form. 

An initial model of a Plummer sphere ((3 = 5 in equa- 
tion 10) is set up with 8000 particles. The satellite is placed 
at the 90% mass radius of the parent. Its orbit is such 
that it is at apocentre, and it has an angular momen- 
tum k times that of the circular orbit with the same en- 
ergy. For k — 0.3, the apocentre (pericentre) of the satel- 
lite orbit is thus approximately 4R (0.1R), where R is the 
scale radius of the parent Plummer sphere. The satellite 
mass is drawn randomly according to equation (11) with 
(/min, /max, rj) = (.05, .2,0). The different values of v and n 
used are summarised in Table 2. The satellite is produced 
by randomly eliminating a fraction (1 — /) of the particles 
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r 

Figure 7. The average profile of the last 7 merger products of 
simulation A (crosses), B (plus signs) and CI (stars) with a fit to 
the inner parts overlaid. The dashed line is a fit to the original 
profile measured in the same way (the fit parameters match a 
Plummer model well). The triangles close to the dotted line are 
a sanity check on the TV-body code (see text). 
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Figure 8. The ratio of T\ji in the parent halo to that of the 
disrupted satellite particles after the merger, as a function of the 
number of mergers, in simulations A (pluses), B (stars) and C[123] 
(dots). 

of the parent and then scaling the positions and velocities 
of the remaining particles. The smoothing length used in 
the A r -body code is a fixed fraction of the simulation box, 
which is not a fixed fraction of the halo scale length. The 
largest value it took in any of the simulations was one third 
of the radius of the innermost point at which the density 
was measured. In simulation A it was at most an eighth of 
this radius. Simulation A was also re-run with the smooth- 
ing length 10 times larger — the result was not significantly 
different. 

After each merger another satellite is generated in the 
same way, and then merged. We measure the profile of the 
halo by starting at the potential centre and binning outwards 



Table 2. The parameters of the TV-body simulations in Sec- 
tion 3.2. The input parameters are (u, k), and the output fitted 
profiles are specified by a. The last column is the value of a pre- 
dicted by equation (8). 



Model 


V 


K 


a 


3*7(1 


A 


0.33 


0.3 


0.66 


0.75 


B 


0.75 


0.3 


1.14 


1.2 


CI 


1.5 


0.3 


1.75 


1.8 


C2 


1.5 


0.6 


1.63 


1.8 


C3 


1.5 


0.05 


1.79 


1.8 



in radius. The 10%, 20%,... radii are determined at the same 
time. Figure 6 shows the evolution of the ratios of these 
decile radii to the half-mass radius r\n. This shows that the 
structure converges after 7 or 8 mergers. Figure 7 shows the 
average profile of the last 7 mergers for models A, B and CI. 
Also shown is the initial Plummer sphere, and a sanity check: 
the profile of the initial model, evolved in isolation for the 
same time as all the mergers added together in simulation 
A. When evolved in isolation the profile changes a little due 
to numerical effects, but much less than when mergers are 
included. The cusp slope a was measured in these models by 
fitting a power law to the profile at r < 0.3r 1( / 2 . The results 
are given in Table 2 and appear as solid lines in Figure 7. 
The general prediction that larger v produces larger a is 
supported, and the actual values of the measured slopes are 
quite close to the prediction of equation (8). 

A comparison of the three models with v = 1.5 (CI, C2, 
C3) shows that the profiles do not differ significantly. Thus 
the final profile is not sensitive to the value of k. 

Cosmological halos are not isolated, as our simulated 
halos are, and they are not in equilibrium in their outer 
parts. Thus we cannot make direct inferences about the 
outer profiles of cosmological halos on the basis of our sim- 
ulations. Isolated merger products are expected to have an 
outer profile with (3 = 4 (Jaffe 1987, Merritt, Tremaine & 
Johnstone 1989). Thus in Figure 7 we see a slow roll-over 
from the inner cusp to an outer envelope with j3 ~ 4, signifi- 
cantly less steep than the initial Plummer model. In cosmo- 
logical halos, f3 ~ 3 is observed near the virial radius for a 
wide range of models. It is not easy to identify the equiva- 
lent of the virial radius in our simulations, but the effective 
(3 of our profiles is about 3 in the range (0.5, 2)riy 2 . The 80% 
radius is at about 2r 1/ / 2 , and the 90% radius (the apocentre 
of the satellite orbits) is at about 2.5r!/ 2 . For comparison, in 
cosmological infall models the turnround radius is at about 
three times the virial radius. 

We also test the 'sinking satellite' assumption of Sec- 
tion 2, which was used in the calculations of Section 3.1. 
If dense satellites really sink to the centre of the parent 
then the particles which originated in the satellite should be 
concentrated towards the centre of the merger product. We 
test this by measuring the half-mass radius of the merger 
product and comparing it with that of the disrupted satel- 
lite particles. The result is shown in Figure 8. The degree of 
concentration of the disrupted satellite particles varies as the 
profile changes. The effect is most pronounced when the ha- 
los have a homogeneous core, and the satellite is much denser 
than the parent in its inner parts. As the cusp develops the 
satellite is disrupted at a progressively larger radius, as we 
argued should happen in Section 2, and the concentration 
settles down at a value of around 2 for all the simulations. 
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4 DISCUSSION 

The analytic argument leading to equation (8) is highly ide- 
alised in its use of a single power law for the cusp profile. 
The semi-analytic model we have used for merging is based 
on the same idea but does not assume a form for the pro- 
file. Nevertheless it is probably not very accurate in detail. 
The iV-body experiment gives us confidence that our sim- 
ple arguments might be close to the truth. The values of 
a, measured from our experiments, agree quite well with 
the prediction of equation (8). Generally the semi-analytic 
experiments give a values higher than predicted, and the 
iV-body experiments give lower values. 

The physical origin we suggest for the cuspy profiles of 



halos h that sufficiently dense satellites can sink intact to 
the halo centre. The details of our merger prescription may 
affect our detailed conclusions, such as equation (8), but the 
robustness of the main physical idea should lead to similar 
behaviour in more realistic models. Figure 8 shows that the 
basic picture of Section 2 is not far wrong, at least for the 
kind of mergers studied here. 

The extent to which this basic picture applies to actual 
cosmological halos remains to be demonstrated. Navarro, 
Frenk and White (1996) introduced the NFW profile for 
CDM halos. Navarro, Frenk and White (1997), in agreement 
with the lower resolution results of Cole and Lacey (1996), 
find that the NFW density profile continues to be a good fit 
to halos from power-law initial fluctuation spectra. The ex- 
pected trend in the mass-density relation of halos is clearly 
visible in their fits: smaller effective v for more negative n. 
According to the arguments of the present work, we would 
expect to see a corresponding trend in profile properties be- 
tween n — and n = —1.5. Nevertheless Navarro, Frenk and 
White (1997) are able to use a single fitting formula for all 
their halos. 

The apparent contradiction may be accounted for in 
a number of ways. Cosmological halos always show a de- 
gree of scatter in their measured properties because of 
substructure — which is absent by design in our experiments. 
More data on the profiles of similar mass halos would give 
better statistical confidence in statements about trends in 
halo properties, in addition, it is difficult to find a consis- 
tent and dynamically meaningful scale for p and r on which 
to base a comparison of halos. In Figure 7 the trend between 
the profile and v is obvious because all the halos line up so 
well in their outer parts. This in turn is a consequence of 
the fact that they are dynamically isolated, and hence the 
half-mass radius and density can be unambiguously defined. 
Since cosmological halos are not isolated, the same compar- 
ison cannot be made in this case. Next, the mass-density 
relation of halos in Navarro, Frenk and White (1997) has 
some scatter in it. If this is a consequence of a real scatter 
in the underlying relation it should produce a scatter in halo 
profile properties in addition to the effects of substructure. 
Finally, our experiments do not properly account for the ex- 
pected systematic dependence on the power spectrum of the 
distributions of satellite mass and orbital parameters. Per- 
haps these could compensate in some way for the difference 
in core structure produced by the different mass-density re- 
lation. Since in our experiments we find no dependence of 
the profile on mass-distribution or orbital eccentricity, and 
since to exactly counter the effects of the mass-density rela- 



tion would require fine tuning, we do not find this explana- 
tion compelling. 

One thing is clear from our experiments though: it takes 
only a few mergers to establish the cusp profile. This means 
that the cusp can adjust to local conditions quite quickly 
and so should not depend in a complicated way on the full 
merger history of a halo. 



REFERENCES 

Aguilar, L., Ostriker, J. P., & Hut, P., 1988, ApJ, 335, 720 

Barnes, J.E., & Hernquist, L., 1991, ApJ, 370, L65 

Cole, S., & Lacey, C, 1996, MNRAS, in press, [astrc-ph/95101 17 
Couchman. H.M .P.. Pearce, F.R., & Thomas, PA., 1996, astro- 



ph/9603116 



Davis, M., & Peebles, P.J.E., 1977, ApJ Supp., 34, 425 

Dehnen, W., 1995, MNRAS, 274, 919 

Dubinski, J., & Carlberg, R.G., 1991, ApJ, 378, 496 

Duncan, M.J., Farouki, R.T., & Shapiro, S.L., 1983, ApJ, 271, 22 

Fillmore, J. A., & Goldreich, P., 1984, ApJ, 281, 1 

Gunn, J.E., & Gott, J.R. Ill, 1972, ApJ, 176, 1 

Hernquist, L., 1990, ApJ, 356, 359 

Hoffman, Y., k, Shaham, J., 1985, ApJ, 297, 16 

Jaffe, W., 1987, Structure and Evolution of Elliptical Galaxies, 

p. 511, ed. de Zeeuw, T., Riedel, Dordecht. 
Kaiser, N., 1986, MNRAS, 222, 323 
McGlynn, T., 1984, ApJ, 281, 13 

Johnstone, D.. 



Merrit, D., Tremaine, S., 
829 

Navarro, J.F., Frenk, C.S. 
270 

Navarro, J.F., Frenk, C.S. 
Navarro, J.F., Frenk, C.S., 
Tormen, G. . Bouchet. F. 
mitted, 



1989, MNRAS, 236, 



& White, S.D.M., 1995, MNRAS, 275, 



& White, S.D.M. 
& White, S.D.M. , 
White, S.D.M. 



istro-ph 9603132 



, 1996, ApJ, 462, 563 
1997, in preparation. 
1996,MNRAS, sub- 



Tremaine, S., etal, 1994, AJ, 107, 634 
van Horner, S., 1957, ApJ, 125, 451 
Villumsen, J.V., 1982, MNRAS, 199, 493 
1994a, AJ, 108, 1398 
1994b, AJ, 108, 1403 
1994c AJ. 108. 1414 



Weinberg, M.D.. 
Weinberg, M.D.. 
Weinberg, M.D.. 
Weinberg, M.D.. 
White, S.D.M. . 
White, S.D.M.. 



1996, 



astro-ph/9607099 



1979, MNRAS, 189, 831 
& Zaritsky, D., 1992, ApJ, 394, 1 
Zhao, H.S., 1996, MNRAS, 278, 488 



This paper has been produced using the Royal Astronomical 
Society/Blackwell Science Tp^X macros. 



© 0000 RAS, MNRAS 000, 0-0 



